Chapter 3A 



N 9 2 - 1 9 l 


Natural Cycles, Gases 


A. R. Douglass, C. H. Jackman 
R. B. Rood, A. C. Aikin, and R. S. Stolarski 
National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, MD 


M. P. McCormick 

National Aeronautics and Space Administration f ’/ 

Langley Research Center [j\ }/ 

Hampton, YA J /flA 

D. W. Fahey r ^ * 

National Oceanic and Atmospheric Administration : L 
Aeronomy Laboratory 
Boulder, CO 





ABSTRACT 


The major gaseous components of the exhaust of stratospheric aircraft are expected to be 
the products of combustion (CO 2 and H 2 O), odd nitrogen (NO, NO 2 , HNO 3 ), and products 
indicating combustion inefficiencies (CO and total unbumed hydrocarbons, THC). The species 
distributions are produced by a balance of photochemical and transport processes. A necessary 
element in evaluating the impact of aircraft exhaust on the lower stratospheric composition is to 
place the aircraft emissions in perspective within the natural cycles of stratospheric species. 
Following are a description of mass transport in the lower stratosphere and a discussion of the 
natural behavior of the major gaseous components of the stratospheric aircraft exhaust. The 
roles of soot and SO 2 are discussed in the next chapter. 


MASS TRANSPORT IN THE LOWER STRATOSPHERE 


The natural cycle of long-lived species is dominated by transport. A qualitative picture 
of the zonal and time averaged stratospheric transport was proposed by Brewer (1) to explain 
the observed sharp decrease in water vapor at the middle latitude tropopause. In this descrip- 
tion, air rises in the tropics and is transported poleward and downward at middle and high lati- 
tudes. Transport from the troposphere to the stratosphere is limited to the tropics, where the 
tropopause temperatures are coldest. The temperatures are cold enough to produce saturation 
with respect to ice, and the ice particles fall. This freeze-dry mechanism thus provides a limit to 
the amount of water vapor that is allowed to enter the stratosphere. Although the water vapor 
measurements are even lower than would be allowed by the average temperatures of the tropical 
tropopause, indicating that the actual mechanism by which air enters the stratosphere is more 
complicated, this qualitative picture of mass transport is consistent with measurements of other 
species. Dobson (2) pointed out that this concept of a circulation (the Brewer-Dobson circula- 
tion) is consistent with ozone observations, which show the highest ozone concentrations in the 
polar stratosphere, away from the tropical source region. Observations of species such as N 2 O, 
which have tropospheric sources and significant loss processes only in the mid to upper strato- 
sphere, are also compatible with the Brewer-Dobson circulation. 

Even though transport processes in the atmosphere are three-dimensional (3-D), two- 
dimensional (2-D) formulations provide a convenient conceptual framework to discuss global 
processes. Furthermore, 2-D models remain at the center of assessment calculations because 
they are computationally manageable. On the most basic level, stratospheric transport can be 
viewed as advection by a large-scale mass transport circulation such as the Brewer-Dobson cir- 
culation discussed in the previous paragraph. Coupled with this advective circulation is large- 
scale quasi-horizontal mixing caused by Rossby waves. The general characteristics of these 
two transport mechanisms are lucidly discussed by Mahlman (3). 

Even though the 2-D framework does offer a useful format in which to discuss global 
transport processes, the mechanics of casting the 3-D circulation into two dimensions are not 
straightforward. Because of this, the best determination of the 2-D advective transport is not 
clear. The simplest idea is to use geometric zonal means about latitude circles to calculate a 
meridional and vertical velocity. This sort of averaging, termed the Eulerian mean, yields a 
meridional flow that is not like the Brewer-Dobson circulation. This difference between the 
expected mass flow and the Eulerian mean arises because of a compensating circulation that is 
driven by the planetary (Rossby) waves. This compensating circulation is a practical manifesta- 
tion of the non-interaction theorem which states that planetary waves have no effect on the mean 
flow in the absence of transience, dissipation, or nonlinearity. While these constraints are never 
strictly met, there is a large compensation between the waves and the Eulerian mean flow, and 
hence, the advective mass circulation is a small residual that must be derived from two larger, 
competing terms. 
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Alternatives to the Eulerian mean have proven very effective in recasting the 2-D prob- 
lem. These are discussed in some detail in chapter 6 of Atmospheric Ozone 1985 (4). 
Theoretically, the most useful framework has proven to be the Lagrangian mean circulation. 
Averages in the Lagrangian mean follow wavy material tubes, and the Lagrangian mean velocity 
tracks this material tube. In the instance of the Eulerian circulation mentioned above, it is found 
that the Lagrangian mean circulation associated with planetary waves is an advective circulation 
directed opposite to the Eulerian zonal mean and with slightly greater magnitude. Dunkerton (5) 
showed that the Brewer-Dobson view of mass transport should be interpreted as a Lagrangian 
mean. Though the concept of the Lagrangian mean is useful theoretically, the necessity of aver- 
aging over wavy material tubes makes a Lagrangian mean model impractical. 

Other estimates of the mean meridional circulation include the transport circulation, the 
diabatic circulation, and the residual circulation. Of these estimates, the residual circulation is 
the most frequently used because it is relatively simple to calculate from transformed Eulerian 
mean equations. These equations are derived by consideration of the large-scale equations of 
motion and combining the heat flux term with the Eulerian advective velocity, to yield a residual 
advective velocity. This procedure effectively captures some, but not all, of the characteristics 
of the noninteraction theorem. Comparisons of the Lagrangian mean and residual circulations 
for simple systems were discussed by Rood and Schoeberl (6), who showed differences 
between the two circulations on the order of 30%. 

Residual circulation models have had considerable success in their representations of 
stratospheric constituent distributions. Calculations for species such as N2O and column O3, 
using both Eulerian mean and residual circulation 2-D models, are in qualitative agreement with 
available global measurements (7). The fact that the residual circulation does offer an approxi- 
mation to the mass flow of the atmosphere is useful in the discussions below. The transport 
effects of the planetary waves, "diffusive transport," constitute a complicated issue that will not 
be considered further here. Both chapter 6 of Atmospheric Ozone 1985 (4) and Andrews et al. 
(8) offer more detailed discussions of transport processes. 

A crucial issue for the High-Speed Research Program (HSRP) is the transport from the 
stratosphere to the troposphere. Exchange processes are discussed in detail in chapter 5 of 
Atmospheric Ozone 1985 (4) and they are considerably more complicated than can be repre- 
sented by 2-D transport descriptions. Accurate representations of convective processes and 
subplanetary scale eddies are essential ingredients of stratosphere/troposphere exchange. Three- 
dimensional models have not effectively dealt with these problems. The data record has only 
begun to identify and address the mechanisms of exchange. Aside from tropical convection, it 
is very likely that significant mass exchange is driven by quasi-horizontal eddies in latitudinally 
and longitudinally confined regions. 

In a time-averaged global domain, the residual circulation can be used to provide at least 
a crude estimate of the mass exchange between the stratosphere and the troposphere. This can 
then be used to estimate stratospheric residence time and, hence, the relative perturbation to the 
background field caused by a pollutant. Consider a stratospheric volume bounded below by the 
tropopause, above by the 50 mb surface, and extending in latitude (4>) from 30°N to 90°N. The 
vertical mass flux, F, is calculated according to 

<(> = 90°N 

F = 2k r e 2 J p 0 w* d(sin<J>) (1) 

4) = 30°N 
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where r e is the radius of the earth (cm), p 0 is the density (g cm- 3 ), and w* is the residual vertical 
velocity (cm s- 1 ) across the tropopause or the 50-mb surface. As noted previously, the estimate 
of mass exchange calculated using the residual circulation is incomplete because eddy transports 
are not included. There is additional uncertainty because the residual circulation is determined 
from the heating rate calculations in the lower stratosphere. The heating and cooling terms are 
nearly balanced, and the net heating from which the residual circulation is derived is small. 
Sensitivity of calculations to small changes in the heating rates is considered by Jackman et al. 
(9). 


The value of F^p is estimated to be 1.4 - 2.4 x 10 17 kg yr- 1 . This estimate is based on 
an evaluation of two residual circulations, RC1 from the NASA/Goddard 2-D model (10) and 
RC2 from the model discussed by Shia et al. (1 1). These may be compared with the estimate of 
Holton (12) for the mass flux through the 100-mb surface of 1.3 x 10 17 kg yr- 1 , and the 
minimum estimate of 1 x 10 17 kg yr- 1 given by Robinson (13). An extreme upper limit of 5.8 x 
1017 kg yr-i is obtained using a circulation RC3, which has larger values of the diabatic heating 
in the lower stratosphere. This is considered an extreme limit because of poor agreement 
between the annual variation of total ozone calculated using the RC3 circulation and Total Ozone 
Mapping Spectrometer (TOMS) measurements (9). 

The downward vertical mass flux across the 50-mb surface is about 9 x 10 16 kg yr 1 for 
RC1, RC2, and RC3. Because the vertical flux out of the layer exceeds the vertical flux into the 
layer, there must be poleward horizontal transport into the region across 30°N. The total mass 
in this region is about 21 x 10 16 kg; the outflow values cited previously imply a turnover time of 
0.9 to 1.5 years for this northern-hemisphere segment of the lower stratosphere. 

This simple model may be used for first-order estimates of the buildup of pollutants in 
the lower stratosphere resulting from aircraft exhaust. For any species with a continuity equa- 
tion may be written for the stratospheric volume described in the previous section 

d%/dt = [Fin Xi - Fout x\ / M + P/M - L% (2) 

where Fj n and F out are the annual mass fluxes of air into and out of the region, is the mass 
mixing ratio of the incoming air-mass flux (Fi n ), M is the air mass of the stratospheric volume 
(kg), P is the total photochemical production rate (kg s- 1 ), and L is the averaged loss frequency 

(s- 1 ). In the presence of aircraft, there is an additional source (P a ircraft)- At steady-state, dyjdl - 
0, and the additional % associated with the aircraft source is 


^aircraft — Faircraft / (Fout + LM) (3) 

An upper limit estimate to x a ircraft can be made by dropping the chemical term (LM). For fuel 
use of 7.7 x 10 10 kg yr- 1 and the outflow estimated above. 


Xaircraft = (7.7 x 10 7 El) / (1.4-2.4 x 10 37 ), (4) 

where El is the emission index (g/kg fuel). The conversion to volume mixing ratios, summa- 
rized in Table 1, is performed by scaling with the ratio of 29 over the mean molecular weight of 
the emitted species. 
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Table 1 . Approximate Perturbations %’ to Stratospheric Trace Gases 
{7.7 x 10 10 kg of fuel/year) 


Species 

El 

(gm/kg fuel) 

Molecular 

weight 

Background 

Perturbation 

co 2 

3100 

44 

350 ppmv 

0.7 - 1.2 ppmv 

GO 

5 

28 

10-50 ppbv 

1.7 - 1.8 ppbv 

H 2 0 

1200 

18 

2-6 ppmv 

0.6 - 1.1 ppmv 

NOy 

15 

46 

2-16 ppbv 

3.0 - 5.2 ppbv 


CARBON DIOXIDE (C0 2 ) 

Carbon dioxide is a product of combustion and one of the major components of aircraft 
exhaust. The El for C0 2 (chapter 1) is roughly 3100 g of C0 2 per kg of fuel. Carbon dioxide 
is chemically inert in the lower stratosphere, i.e., there are no significant photochemical produc- 
tion or loss processes. Measurements of the mixing ratio profile indicate little height depen- 
dence (14,15), which is expected in view of the steadily increasing tropospheric values. 
Current tropospheric measurements indicate a volume mixing ratio of about 350 ppmv (16). 
Assuming that this value is representative of the lower stratosphere, there are approximately 1 x 
IOI 4 kg in the bounded area described in the previous section. For an estimated fleet fuel usage 
of 7.7 x 10 10 kg/yr, this corresponds to the source of C0 2 of 2 x 10 11 kg yr-k The average 
increase in C0 2 over this volume is 0.7 to 1.2 ppmv (Table 1), less than 1% of the background 
level. 


CARBON MONOXIDE (CO) 

Carbon monoxide reacts in the lower stratosphere with the hydroxyl radical. The photo- 
chemical lifetime, taken from the NASA/Goddard 2-D model (10), is given in Figure 1. 
Although short lived in the upper stratosphere, the lifetime of CO in the lower strato-sphere (1 
to 4 years) is comparable to its turnover time (1 to 3 years). Carbon monoxide provides a mea- 
sure of combustion inefficiency, and the estimated El values are all less than 5 g of CO per kg 
of fuel. Although this source of CO (< 3.8 x 10 8 kg/yr for 7.7 x 10io kg of fuel) may be compa- 
rable to the natural source of CO (about 6.5 x 10 8 kg/yr for the bounded region, largely a result 
of methane oxidation), the long lifetime indicates that the perturbation will be controlled by mass 
flow through the region. If chemical effects are neglected, the expected perturbation to the 
volume mixing ratio is less than 2 ppbv, compared with the natural level in the lower strato- 
sphere of 10 to 50 ppbv (Table 1). 


WATER VAPOR (H 2 0) 

Water vapor is a product of combustion and an important component of aircraft exhaust. 
The normal H 2 0 mixing ratios in the lower stratosphere are low (3-5 ppmv), as will be dis- 
cussed further in this section. The perturbation is estimated to be 0.6 - 1.1 ppmv, which is 
significant compared with the background values (Table 1). Because of its importance to the 
lower stratosphere, a more thorough examination of available data is warranted. Water vapor is 
the major source for odd hydrogen radicals in the stratosphere, generated through reactions with 
excited atomic oxygen O(tD). Odd hydrogen reactions contribute to the ozone loss throughout 
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the stratosphere and are most frequent in the lower stratosphere and mesosphere (17). They 
also modify the ozone destruction efficiency of odd nitrogen and odd chlorine compounds, 
making chlorine more effective and nitrogen less effective. The amount of water vapor incopjo- 
rated in aerosols is potentially critically important for the efficiency of heterogeneous reactions 
(18,19,20), which may directly affect the balance of odd nitrogen and odd chlorine species in 
the lower stratosphere. The formation of polar stratospheric clouds, either nitric acid trihydrate 
or ice clouds, depends on the stratospheric temperature and on the amount of water vapor avail- 
able. The most important photochemical process that affects stratospheric water is oxidation of 
methane (CH 4 ), which produces about half of the water vapor in the upper stratosphere. A 
detailed discussion of the photochemistry of CH 4 , the production of stratospheric H 2 O, and the 
role of molecular hydrogen (H 2 ) is provided by LeTexier et al. (21). 
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Figure 1. The photochemical lifetime of CO (years) (March) taken from the NASA/GSFC 2-D model. 


Water vapor measurements have long been considered a useful indicator of stratospheric 
transport. Early sonde measurements of the water vapor profile showed a rapid decrease in the 
mixing ratio profile above the midlatitude tropopause. The measured values were substantially 
lower than the tropopause saturation mixing ratio at middle latitudes. The Brewer-Dobson cir- 
culation was proposed as an explanation for these observations (1,2). Rising air enters the 
stratosphere only in the tropics, where there is rising motion, and is transported poleward and 
downward at middle and high latitudes. The stratospheric water vapor mixing ratio is limited by 
the cold temperatures of the tropical tropopause. More recent measurements (22) have demon- 
strated the existence of a hygropause, a minimum in water vapor mixing ratio above the 
tropopause at 19-22 km between 30°N and 30°S. 
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Various techniques have been applied to the measurement of stratospheric water. These 
were reviewed by the WMO in 1986 (4). Although there are measurements of individual pro- 
files that use assorted techniques, until recently balloon-borne frost point hygrometer measure- 
ments provided the most complete picture of annual behavior. In Figure 2, taken from 
Mastenbrook and Oltmans (23), the annual variation of water vapor (mass mixing ratio) is given 
for 1964-76 at Washington, DC (39°N, 77°W). The annual cycle near 50-100 mb is much 
weaker compared with the annual cycle near 150 mb. Similar behavior is seen for a data set 
taken at Boulder, CO (40°N, 105°W) for 1981-89 (Figure 3). In the lower stratosphere, near 
70 mb, there is still a seasonal cycle, with a maximum value in the second half of the year. The 
seasonal motion of the tropopause is clearly visible, reaching a maximum altitude (i.e., highest 
H 2 O) in August to September. The increase in mixing ratio with altitude indicates the produc- 
tion of H 2 O from CH 4 . 

The first global measurements of water vapor were provided by the Limb Infrared 
Monitor of the Stratosphere (LIMS) experiment for the period November 1978 through May 
1979 (24). The vertical resolution of the measurements is 4 km. At 50 mb, there is a pro- 
nounced latitudinal gradient, as shown in Figure 4, consistent with upward motion in the trop- 
ics. The gradient is much weaker at 10 mb. The magnitude of the annual variation is consistent 
with that derived from the frost point hygrometer measurements shown in Figures 2 and 3. 
Zonal mean profiles for each month during the LIMS period are given in Figure 5(a-g). These 
measurements are in general agreement with rising motion in the tropics, downward motion at 
the high latitudes, and a high- altitude source. 



Figure 2. Annual variation of stratospheric water vapor (mass mixing ratio) based on balloon-borne 
frost point hygrometer measurements at Washington, DC, for the period 1974-76 (Mastenbrook and 
Oltmans [23]). 
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Figure 4. LIMS water vapor latitude versus time cross section at 50 mb (0.25-ppmv intervals, 
5-day means). 
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Figure 5. LIMS monthly zonal mean water vapor pressure versus latitude cross sections for 
(a) November and (b) December 1978 and (c) January, (d) February, (e) March, (f) April, and 
(g) May 1979 (0.5-ppmv contour intervals). 
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The Stratospheric Aerosol and Gas Experiment (SAGE II) aboard the Earth Radiation 
Budget Satellite was launched in 1984. The vertical resolution of this solar occultation instru- 
ment is 1 km. The SAGE II measurements provide a much longer data record than previously 
available, and furnish a picture of seasonal and year-to-year variability in the upper troposphere 
and stratosphere over much of the globe. The water vapor data, validation studies, and initial 
results are discussed by Rind et al. (25,26), Chiou et al. (27), Chu et al. (28), Larsen et al. 
(29), McCormick et al. (30), and Oltmans et al. (31). Validation studies conducted after the 
Sage II launch have included comparisons with correlative balloon-borne frost point hygrometer 
and a limited number of aircraft Lyman-alpha flights. Figures 6(a) and 6(b) are examples of 
comparisons with frost point measurements over Boulder and Lyman-alpha measurements over 
southern high latitudes, respectively. 
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Figure ©(a). Comparison of SAGE li measurements with frost point hygrometer measurements near 
Boulder CO. 
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Figure 6(b). Comparison of SAGE II measurements with aircraft Lyman-alpha measurements at 
southern middie latitudes. 


A time series of the average water- vapor-mixing ratio between 30° N and 40° N is given 
in Figure 7. There is excellent agreement between the seasonal pattern in the SAGE II data and 
the 7 years of balloon measurements made over Boulder (Figure 3). Both show a minimum in 
early winter near the altitude range 14-15 km. The minimum has its lowest value in March, and 
thereafter weakens and shifts to higher altitudes. The balloon data show that the magnitude of 
the seasonal variation near 100 mb at mid to high latitudes is smaller than might be expected 
based on the LIMS data, but the phases of the variations are similar. 

The altitude-latitude cross section of SAGE II water vapor mixing ratios are given for 
May 1987 in Figure 8. The SAGE measurements are also consistent with the Brewer-Dobson 
view of mass transport and a high altitude source. These measurements compare well with the 
measurements in Figure 5(g) from the LIMS instrument. 
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Figure 8. The latitude-altitude distribution of SAGE II measurements of water vapor (volume 1 mixing 
ratio, ppmv) for May 1987. 
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ODD NITROGEN 


Odd nitrogen species (primarily NO and NO 2 ) constitute an important trace component 
of aircraft exhaust. The anticipated perturbation to total odd nitrogen is significant compared 
with the normal background values, to be discussed later. Odd nitrogen species (total odd 
nitrogen = NO y = N+N 0 +N 0 2 +N 0 3 + 2 *N 2 05 +HN 0 3 + HN0 4 +C10N0 2 ; NO x = 
NO+NO 2 .+NO 3 ) are extremely important to the stratospheric ozone balance. Odd nitrogen is 
computed to account for about 50% of the total ozone loss below 40 km and between 60 and 
90% of the ozone loss in the lower stratosphere (17). Although there are many measurements 
of individual odd nitrogen species, there are few measurements that can be interpreted as pre- 
senting a global picture of total odd nitrogen. The natural cycles must be derived from a combi- 
nation of models and data. Following is a discussion of sources of stratospheric odd nitrogen, 
including nitrous oxide oxidation, galactic cosmic rays (GCRs), lightning (with upflux from the 
troposphere to the stratosphere), solar proton events (SPEs), relativistic electron precipitations 
(REPs), and dissociation of N 2 in the thermosphere (with downflux through the mesosphere to 
the stratosphere). Nuclear explosions, specifically above-ground detonations with energies 
above 1 megaton, also can produce large amounts of odd nitrogen in the stratosphere. The only 
loss processes for total odd nitrogen are recombination of N and NO in the upper stratosphere 
and mesosphere, and transport to the troposphere. 

Natural Sources of Odd Nitrogen 

The major source of stratospheric odd nitrogen results from nitrous oxide oxidation 
(N 2 0 + OOD] ->NO + NO); several groups have computed globally averaged annual production 
rates (32-36). An example of the model computation of the odd nitrogen production rate from 
the oxidation of N 2 0 for March (dashed line) is given in Figure 9. This production rate was 
taken from a 2-D model computation wherein the N 2 Q levels were constrained by stratospheric 
and mesospheric sounder (SAMS) measurement, and the O 3 levels were constrained by mea- 
surement from the solar backscatter ultraviolet (SBUV) data [see Figure 9a, Jackman et al. 
(36)]. The global annually averaged production of odd nitrogen from this source is about 2.6 x 
10 34 molecules yr 1 and is compared with other odd nitrogen sources in Table 2. 


Table 2. Comparison of Odd Nitrogen Sources and Sinks in the Stratosphere 


Source 


Nitrous oxide oxidation 
N 2 0 + 0(!D) - > NO = NO 

26. 

Galactic cosmic rays 

3.7 (solar minimum) 

2.7 (solar maximum) 

Lightning 

11 . 

Solar proton events 

1.5 (maximum in 1972) 

Thermospheric downflux 

7 

Relativistic electron precipitations 

? 

Nuclear explosions 

24. (1961 & 1962 nuclear tests) 

High speed civil transport planes 

SINK 

14. (emission index =15) 

Reforming of molecular nitrogen 
N + NO - > N 2 + O 

8.4 

Transport to the troposphere 

32. 
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Galactic cosmic rays produce odd nitrogen in the lower stratosphere and upper tropo- 
sphere. The earth's magnetic field influences the GCRs to a certain extent, such that more 
GCRs affect the atmosphere at higher geomagnetic latitudes. Galactic cosmic rays produce odd 
nitrogen through dissociation or dissociative ionization processes in which N 2 is converted to 
N( 4 S), N( 2 D), or N+ (33,37,38). Rapid chemistry is initiated after N 2 dissociation, and most 
of the atomic nitrogen is rapidly converted to NO and NO 2 . The GCR mean production rate of 
odd nitrogen [solid line, taken from Figure 13 of Jackman et al. (36)] is compared with the oxi- 
dation of N 2 O source in Figure 9. A solar cycle variation is apparent in the GCR flux. The 
stratospheric component of the odd nitrogen production is estimated to be a maximum of 3.7 x 
10 33 molecules yr-i during solar minimum and a minimum of 2.7 x 10 33 molecules yr 1 during 
solar maximum (33). 

Natural production of odd nitrogen is dominated in the middle and upper stratosphere by 
N 2 O oxidation and in the lower stratosphere at the higher latitudes by the GCRs. Since the odd 
nitrogen family has a lifetime of months in the middle and lower stratosphere, transport of odd 
nitrogen created at higher altitudes and lower latitudes is significant, and thus the GCR source 
of odd nitrogen has been computed to increase odd nitrogen in the lower stratosphere at high 
latitudes by only about 10% (37). 



-90S -60 -30 0 30 60 


LATITUDE (deg) 


Figure 9. Odd nitrogen production (number cm 3 sec' 1 ) caused by GCRs (solid line, from Figure 13 of 
Jackman et al. (36) and oxidation of nitrous oxide (dashed line, from Figure 9 of Jackman et al" (36). 
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Lightning is responsible for the production of significant amounts of odd nitrogen (of 
order 3 x 10 35 molecules yr- 1 ) from N 2 in the free troposphere (35, 39-43). This source, along 
with others of tropospheric NO y , can reach the stratosphere through the transport of tropo- 
spheric air across the tropopause. For example, the model source of lightning designated "H2" 
by Ko et al. (35) transports 1.1 x 10 34 molecules of NOy yr- 1 to the stratosphere when included 
in the NASA/Goddard 2-D model. Lower stratospheric background mixing ratios are increased 
substantially by including this lightning source. However, preliminary analysis of the STEP 
(Stratosphere Troposphere Exchange Project) data near the tropical tropopause show exceed- 
ingly low values of odd nitrogen (<0.5 ppbv). The correlation of odd nitrogen with N 2 0 found 
on the aircraft missions to both the Antarctic (AAOE, Airborne Antarctic Ozone Experiment) and 
the Arctic (AASE, Airborne Arctic Stratospheric Expedition) is also consistent with the view 
that air entering the stratosphere contains low concentrations of odd nitrogen (44a,44b). 

The production of odd nitrogen species by SPEs was predicted and shown to affect 
upper- stratospheric ozone (45). Since solar protons generally are at lower energies than GCRs, 
they produce NO at higher altitudes; they are constrained more significantly by the earth's mag- 
netic field and affect the atmosphere only in the polar cap region (geomagnetic latitudes greater 
than about 60°). SPEs are sporadic, with durations up to several days; more SPEs occur close 
to the solar maximum. 

The increase in polar NO after the July 1982 SPE was inferred to be about 6 x 10 14 NO 
molecules cm- 2 (4.1 x 10 32 molecules NO y ) at polar latitudes from the SBUV instrument (46), 
in good agreement with the calculated NO increase of 7 x 10 14 NO molecules cm- 2 (4.8 x 10 32 
molecules NO y ) in the polar cap (47). The August 1972 SPE was one of the biggest events in 
the past 30 years, and computations revealed that a substantial production of odd nitrogen (1.5 x 
IQ 33 molecules NO y ) was associated with this event at polar latitudes in the middle to upper 
stratosphere (47). 

The influence of extreme ultraviolet, auroral electrons, and photoelectrons on the odd 
nitrogen budget of the stratosphere, through transport of odd nitrogen from the thermosphere 
downwards through the polar mesosphere to the upper stratosphere, has been studied by several 
groups (33,38,48-56). EUV, auroral electrons, and photoelectrons are capable of dissociating 
N 2 to form atomic nitrogen in the thermosphere. Transport of this odd nitrogen to the meso- 
sphere and upper stratosphere is possible, but certain conditions must be present. The lifetime 
of odd nitrogen in the sunlit mesosphere is short, and it is only during the long period of polar 
night at high latitudes, when a period of several weeks of darkness is typical, that significant 
downward transport is possible. Solomon et al. (51) undertook a detailed 2-D model study of 
thermosphere-middle atmosphere coupling. Enhancements in the odd nitrogen distribution in 
the mesosphere and upper stratosphere were found to occur when the thermospheric production 
of odd nitrogen was included [compare Figures 8 and 17 of Solomon et al. (51)]. These 
enhancements of NO x in the middle atmosphere, caused by thermospheric downflux of odd 
nitrogen, are especially significant in the hemisphere that was most recently polar night. 
Measurements by the LIMS of NO 2 have also indicated that odd nitrogen may be enhanced in 
the mesosphere and, to a lesser degree, in the upper stratosphere during polar night (55). The 
magnitude of the thermospheric downflux of NO x into the stratosphere is highly uncertain and 
variable, but probably small. 

In the past 15 years, it has been proposed that REPs make important contributions to the 
polar odd nitrogen budget of the mesosphere and upper stratosphere (57-62). The frequency 
and flux spectra of these REPs is still under discussion. Baker et al. (59) show evidence of 
large fluxes of relativistic electrons at geostationary orbit measured by the Spectrometer for 
Energetic Electrons (SEE) instrument on board spacecraft 1979-053 and 1982-019. REPs, 
which are actually depositing energy into the middle atmosphere, have been measured by 
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instruments aboard sounding rockets (63). These rocket measurements have typically indicated 
much smaller fluxes of relativistic electrons than measured by the SEE instrument. REPs that 
display the large fluxes measured by Baker et al. (59) are computed to show an important influ- 
ence on the NOy budget in the lower mesospheric-upper stratospheric region by Callis et al. 
(62). However, REPs with the smaller fluxes measured by Goldberg et al. (63) would cause a 
relatively insignificant change in odd nitrogen amounts in the middle atmosphere at this produc- 
tion rate. More work is necessary to determine which REP events are more typical. The global 
magnitude of the REP source of stratospheric odd nitrogen is also highly uncertain but probably 
small. 

Sources of Odd Nitrogen Related to Human Activity 

The natural source of odd nitrogen is increasing because the source gas, N 2 O, is 
increasing. The global rate of increase of N 2 O is about 0.7 ppbv yr-l; the global average 
concentration is currently above 308 ppbv (16). At this rate of increase, assuming a linear 
increase in N 2 O at all altitudes, the natural source of odd nitrogen will increase by 10% in 40 
years. For comparison, the source of NOy from stratospheric aircraft for with El of 15 and 7.7 
x 10 10 kg of fuel yr - 1 is 1.4 x 10 34 molecules yr 1 , nearly 40% of the natural source. 

Atmospheric nuclear explosions, affecting the 7-to 30-km regime, cause sufficiently 
high temperatures in the fireball that thermal decomposition of N 2 and O 2 occurs, followed by 
the formation of NO (64-66). The nuclear tests in 1961 and 1962 by the Soviet Union and the 
United States were thought to be responsible for the production of 2.4 x 10 34 molecules of odd 
nitrogen (16). 

Sinks of Odd Nitrogen 

There are two major sinks for odd nitrogen in the stratosphere: the recombination of N 
and NO in the middle to upper stratosphere, to form N 2 and O, and transport to the troposphere. 
Nitric oxide (NO) is photodissociated easily in the upper stratosphere, forming atomic nitrogen 
(N) and atomic oxygen (O). Some of this N then reacts with NO, to form molecular nitrogen 
(N 2 ); the rest reacts with O 2 to reform NO. The lower stratospheric odd nitrogen abundance is 
principally regulated by transport of odd nitrogen to the troposphere, with subsequent rainout of 
the HNO 3 andH 02 N 02 constituents [see Jackman et al. (36)]. 

The global annually averaged loss for the first process is about 8.4 x 10 33 molecules yr- 
1 (Jackman and Douglass, 2-D model computation, 1990). The remainder of the globally aver- 
aged production is balanced by the transport of stratospheric odd nitrogen to the troposphere 
(4.0 x IO 34 - 0.8 x 10 34 = 3.2 x 10 34 ), the dominant loss mechanism for stratospheric odd 
nitrogen. Thus the amount of odd nitrogen in the lower stratosphere depends critically on the 
circulation in the lower stratosphere and stratosphere/troposphere exchange. 

Odd Nitrogen Distribution 

The global distribution of total odd nitrogen may be calculated using a 2-D model, 
keeping in mind that the model transport in the upper troposphere/lower stratosphere is impor- 
tant to determining distribution and residence time of total odd nitrogen. Although there are no 
global measurements of total odd nitrogen, Callis et al. (67) estimate total odd nitrogen from 
nighttime LIMS measurements of NO 2 and HNO 3 . Using a model, they show that this sum 
represents 68 , 78, 85, and 90% of total odd nitrogen at, respectively, 28, 34, 37, and 40 km. 
The estimate from LIMS for March (Figure 10) can be compared to the calculations of the 
Jackman and Douglass model (Figure 1 1), or to several model distributions given by Jackman 
et al. (7). Although the data represent a lower limit because all species are not included, notably 
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N 2 O 5 , the model values are somewhat lower than the measurements in the lower stratosphere. 
It should be noted that the lower stratospheric LIMS values of NO 2 depend on the NO 2 clima- 
tology used in the retrievals resulting from a loss of NO 2 signal from altitudes below 40 mb 
(68). Thus the LIMS data can provide constraints on the odd nitrogen species but do not pro- 
vide a direct measure of NOy. 

In situ measurements of total odd nitrogen taken aboard the NASA ER-2 aircraft provide 
several points for direct comparison between models and data. The ambient total odd nitrogen 
is measured by catalytic conversion of the higher-oxide species to NO and subsequent detection 
of NO in the chemiluminescent reaction with reagent O 3 . All odd nitrogen species of impor- 
tance in the stratosphere are included in this technique (69). Although measured odd nitrogen in 
the lower stratosphere and upper troposphere has substantial natural variability on a wide range 
of scales, the data are sufficient to derive representative averages at some locations and seasons. 
Extensive measurements in the winter hemispheres near 55°S and 60°N show odd nitrogen 
values less than 0.5 ppbv in the troposphere, up to about 10 km altitude, increasing to about 3 
ppbv near 16 km and to about 8 ppbv near 20 km (70,71). Values of 3 to 8 ppbv between 19 
and 21 km were measured near 40°N in October (72). Measurements near 10°S in January and 
February gave tropospheric odd nitrogen values similar to those measured at high latitudes but 
significantly lower values in the lower stratosphere. Odd nitrogen values near 10°S increase 
from about 0.5 ppbv near the tropopause to about 2.5 ppbv near 21 km (73). The model values 
in Figure 1 1 underestimate the measured total odd nitrogen. Previous ER-2 NO measurements 
have suggested this discrepancy (4,74). The ATMOS (Atmospheric Trace Molecule 
Spectroscopy) measurements of total odd nitrogen are as much as 40% larger than the results 
from some model calculations for the lower stratosphere (75). 
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Figure 10. LIMS nighttime HN0 3 + N0 2 for March 1979, a lower limit for total NO y . 


50 





00 


50 E 

£s£ 

LLI 

40 Q 

D 

t 

30 b 

< 


20 

10 

0 


x 

o 

ir 

o. 

a. 

< 


Figure 11. A 2-D model calculation of total odd nitrogen (March). 


Partitioning of Odd Nitrogen Species 

The lifetimes for interchange reactions among key odd nitrogen species, calculated from 
the NASA/Goddard 2-D model, are given for March equinox in Figure 12. Once HNO 3 is 
formed in the lower stratosphere, it is relatively long lived (> 50 days for middle and high lati- 
tudes). The lifetime of N 2 O 5 is short, except in high latitudes in winter. The equilibrium 
between NO and NO 2 is established quickly; thus, the fractions of NO and NO 2 in the aircraft 
exhaust are not important for assessment calculations. Some studies have likewise shown an 
insensitivity to the relative partitioning between NO x and HNO 3 in stratospheric injections, but 
this is not likely to be the case for tropospheric emissions. 

However, the current understanding of odd nitrogen chemistry in the lower stratosphere 
does not explain all of the observations. For example, measured values of wintertime HNO 3 
exceed model calculated values (35). Some investigators have suggested that the heterogeneous 
reaction of N 2 O 5 on polar stratospheric clouds (PSCs) in the polar vortex may produce the 
observed values (76). Rood et al. (77) have found, with 3-D model calculations, that there is 
insufficient mixing between the polar vortex and middle latitudes to account for the HNO 3 
deficit. Measurements of reaction rates show that N 2 O 5 + H 2 0 reactions occur on sulfuric acid, 
even though some of the other processes that occur on PSC surfaces (e.g., HC 1 + CIONO 2 ) do 
not (18). Hofmann and Solomon (78) suggested that such heterogeneous reactions on back- 
ground aerosols could account for the observations. These processes are important for under- 
standing the natural balance of HNO 3 and NO x . The potential impact of aircraft exhaust on 
stratospheric ozone is sensitive to this balance. 
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Figure 12. Lifetimes for interconversions among odd nitrogen species (March). 


HYDROCARBONS 

The lower troposphere contains a wide variety of hydrocarbons. Most of these are 
destroyed by reactions with the hydroxyl radical. Only those species with photochemical life- 
times longer than a few months reach the tropical tropopause in sufficient concentration to affect 
stratospheric chemistry. This restricts the population to methane (CH 4 ) and the nonmethane 
hydrocarbons (NMHC) acetylene (C 2 H 2 ), ethane (C 2 H 6 ), and propane (QjHg). 

The methane distribution is fairly well established on several spatial scales. Data from 
the period 1983-87, from the surface and mid-troposphere, indicate that CH 4 is increasing by 
about 1% per year (16). During the 1980s, the average free tropospheric mixing ratio in the 
northern hemisphere was about 1.7 ppm, while in the southern hemisphere it was somewhat 
less, 1.6 ppm. The stratospheric distribution has been determined by several groups measuring 
primarily at midlatitudes and using balloon-borne instruments. The SAMS measurements pro- 
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vide a global picture of the behavior of CH 4 above 20 mb. The maximum stratospheric mixing 
ratio is seen in the tropical lower stratosphere; the mixing ratio decreases poleward and with 
altitude. A detailed description of the data and comparison with model calculations are given by 
Jones and Pyle (79). 

Nonmethane hydrocarbons have considerable seasonal and spatial variation in the tropo- 
sphere and we may expect similarly large variations in the lower stratosphere as well. Average 
lifetimes for C 3 H 8 , C 2 H 2 , and C 2 H 6 are 0.06 yr, 0.1 yr, and 0.25 yr, respectively, but vary 
according to locale and may be much shorter in the tropics. Based on limited observations, we 
may describe the distribution of NMHCs as follows. The NMHCs increase from the South 
Pole towards the equator, (80-84), peaking at northern midlatitudes. Typical tropospheric 
northern midlatitude values for the NMHCs are 0.2 ppb acetylene, 2 ppb ethane, and 0.6 ppb 
propane. The seasonal cycle is dominated by variations in OH. In the Arctic, ethane has a 
maximum mixing ratio in the winter; and is 2 times less in summer. The annual variation of 
propane is greater (85). In Antarctica ( 86 ), ethane has a winter value of 0.45 ppb, with a 
summer mixing ratio of 0.3 ppb. The corresponding figures for acetylene are 20 and 10 ppt, 
respectively, and the propane content varies from 0.05 to 0.1. ppb. 

Vertical distributions of the NMHCs in the stratosphere have only been measured at 
northern midlatitudes. With the exception of a limited number of satellite measurements, infor- 
mation has been by balloon-borne samplers (83,87-90) or by analysis of solar absorption 
spectra taken from the ground or aircraft (91-93). The vertical profile of ethane is nearly con- 
stant to 10-km altitude, with a mixing ratio between 1.5 and 2 ppb. The mixing ratio drops by 
an order of magnitude at 15 km and by 2 orders of magnitude at 25 km. This vertical distribu- 
tion is confirmed by analysis of spectra profiles of solar infrared radiation (94). The measured 
vertical profiles of propane are confined to altitudes below 20 km. At 10 km, mixing ratios are 
between 0.2 and 0.8 ppb; at 17 km, the mixing ratio is 0.1 ppb. Acetylene, as measured at 10 
km by direct sampling, has a value between 0.1 and 0.23 ppb. By 16 km, the mixing ratio is 
0.02 ppb. The vertical distribution of molecules like ethane could be important in the chemistry 
of the stratosphere by serving as a tracer for the OH and Cl radical profiles (95). The break- 
down of the NMHCs can modify the photochemistry of the lower stratosphere and upper tropo- 
sphere (96-99). 

The lifetimes of NMHCs can be significantly shorter in the stratosphere than in the tro- 
posphere as a result of reactions with Cl and are often less than 1 month. If this estimate of the 
lifetime is correct, the magnitude of the increase in the concentrations of the hydrocarbons in 
aircraft exhaust is controlled by photochemical loss rather than transport. Pollutant molecules 
are destroyed by reactions with OH and Cl, the perturbation is limited directly by the amount of 
pollutant in the aircraft exhaust. An important role of CH 4 and NMHCs in the lower strato- 
sphere is to convert Cl to HC1 and to produce ozone by reactions of NO with oxidation prod- 
ucts, e.g., HO 2 and CH 3 O 2 . The rate of oxidation of background CH 4 in the lower strato- 
sphere (1-5 x 10 10 molecules cm - 3 s- 1 ) is much larger than the source of hydrocarbons in the air- 
craft exhaust (about 1 x 10 8 molecules cm - 3 s 1 if mixed in the northern hemisphere over an 
altitude thickness of 5 km). The picture could be somewhat more complicated by addition of 
reactions to form nitrates such as peroxy acetyl nitrate (PAN) in the high NO x environment of 
the aircraft plume. However, although PAN is somewhat longer lived than the hydrocarbons, 
the summer lifetime is less than 3 months, and the maximum perturbation is still controlled by 
photochemistry. 
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CONCLUSION 


The species considered here, most obviously the odd nitrogen species, influence the 
photochemical balance of the lower stratosphere. Perturbations to these species may have a 
direct impact on the ozone concentration in the lower stratosphere, even though time scales for 
photochemical processes such as ozone destruction are in the range of several years. For the 
species that are long lived in the lower stratosphere, the magnitudes of the possible perturba- 
tions to the background values of the components of stratospheric aircraft exhaust are controlled 
by the mass flow through the region in which the aircraft exhaust is mixed. This underscores 
the importance of understanding stratosphere/troposphere exchange and the circulation in the 
lower stratosphere where, it has been proposed, the high-speed civil transport will fly. 
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